Probing high-energy electronic excitations using inelastic neutron scattering 
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High-energy, local multiplet excitations of the d-electrons are revealed in our inelastic neutron 
scattering measurements on the prototype magnetic insulator NiO. These become allowed by the 
presence of both non-zero crystal field and spin-orbit coupling. The observed excitations are con- 
sistent with optical, x-ray, and EELS measurements of d-d excitations. This experiment serves as 
a proof of principle that high-energy neutron spectroscopy is a reliable and useful technique for 
probing electronic excitations in systems with significant crystal field and spin-orbit interactions. 

PACS numbers: 78.70.Nx, 71.70.Ch, 71.70.Ej 



Electrons in strongly correlated systems are character- 
ized by the dual nature of their spatial wave functions; 
sometimes they are considered localized, but other times 
delocalized. This dichotomy makes it difficult for tradi- 
tional band theory to describe physical properties of these 
scientifically interesting and technologically relevant ma- 
terials. Often a model based on local interactions, such 
as a Hubbard Hamiltonian, is used to capture the main 
physics of these systems, but distilling the numerous in- 
teractions present in a condensed matter system into a 
simple model requires information from experiments. For 
this reason, various types of electron spectroscopy tools 
have played significant roles in elucidating various inter- 
actions such as crystal field (CF) and spin-orbit (SO) 
coupling. These include optical, x-ray absorption, pho- 
toemission, and electron energy loss spectroscopy (EELS) 
as well as inelastic x-ray scattering (IXS) in recent years. 
Here we show that inelastic neutron scattering (INS) can 
be added to this arsenal of spectroscopic tools for probing 
electronic structure of correlated electron systems. 

Neutron spectroscopy has been used to study crystal 
field splittings in rare-earth insulators [l] and metals [2- 
With the availability of copious amounts of epither- 
mal neutrons from spallation neutron sources, even inter- 
multiplet transitions in f-electron systems have been ob- 
served with INS (gI, |7|. Despite earlier theoretical pre- 
dictions of using INS in studying interband transitions 
in semiconductors [8], no experimental studies of elec- 
tronic excitations above 1 eV have been reported to date. 
Only recently, it was shown that quantitative informa- 
tion above 1 eV can be obtained from direct geome- 
try spectrometers [9]. We note that novel IXS methods 
have been introduced utilizing either quadrupole tran- 
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such charge-sensitive probes are not ideal for investigat- 
ing doped Mott insulators, such as cuprate superconduc- 
tors, because charge carrier contributions dominate the 
excitation spectrum in the energy range of interest. On 
the other hand, charge-neutral neutrons do not couple to 
charge degrees of freedom, and could be very valuable for 
measuring d-d excitations in such metallic compounds. 

In this Letter, as a first step towards this goal, we re- 
port our observation of localized d-d excitations in NiO 
using INS. NiO is a prototypical Mott insulator with a 
large charge-transfer gap of about 4 eV. This material has 
been studied extensively using various spectroscopic tech- 
niques, and its electronic structure is quite well known. 
In particular, many novel spectroscopic tools have been 
recently applied to study NiO d-d excitations as a test 
case, since the NiO d-d excitations are well character- 
ized [l5|. For example, spin-polarized EELS [IgI, soft 
x-ray resonant inelastic x-ray scattering (RIXS) [12I. [l3|. 
hard x-ray RIXS [MK non-resonant IXS [10], and angle- 
resolved EELS [l7| have all been used to detect d-d exci- 
tations in NiO. The benefit of studying NiO is that there 
is no ambiguity with regard to the identification of d-d 
excitations, thus allowing us to focus on the mechanism 
of the spectroscopy. 

In our INS investigation of NiO, d-d excitations at 1.0 
eV and 1.6 eV have been observed, which is the high- 
est energy solid-state excitation observed with INS to 
date f?\ . These transitions occur via magnetic dipole and 
higher order scattering operators, respectively, giving 
strong momentum dependence to their spectral weights. 
We find that either CF splitting or SO coupling is nec- 
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FIG. 1: (Color online) Panel a) is the magnetic intensity map 
of energy- momentum space obtained on MAPS using incident 
neutron energy of 3 eV at room temperature. The background 
intensity due to phonons has been subtracted as discussed in 
the text, b) illustrates the theoretical calculation discussed in 
the text. Panels c) and d) illustrate integrated cuts in momen- 
tum and energy as noted. The open circles are experimental 
data and the solid (red) lines are theoretical calculations. As 
described in the text, the overall amplitude is the only ad- 
justable parameter in the calculation. 



essary for the observed transitions, suggesting that one 
could determine the strengths of these interactions using 
momentum dependent INS spectrum. We also studied 
single magnon density of states, which peak up near the 
zone boundary energy. In contrast to recent soft x-ray 
RIXS studies, we do not find atomic spin flip scattering 
at the expected energy range [18]. Our observation thus 
opens up an exciting new field of electron spectroscopy 
using INS at quite an opportune time; recently commis- 
sioned Spallation Neutron Source will undoubtedly pro- 
vide an excellent venue for such studies. 

Total 250 g powder sample of NiO (Alfa Aesar 99 %) 
was studied on the MAPS and MARI direct geometry 
time of flight spectrometers located at ISIS (Rutherford 
Appleton Labs, U.K.). On MAPS, an incident energy of 3 
eV was selected by spinning the nimonic chopper at 100 
Hz and the "A" Fermi chopper at 600 Hz. On MARI, 
an incident energy of 500 meV was used with a 50 Hz 
nimonic and a 600 Hz "A" chopper. All measurements 
were done at room temperature. Note that the Neel tem- 
perature of NiO is well above room temperature. 

In Fig. 1(a), the intensity map of the excitation spec- 
trum of NiO obtained with Ei=3 eV is plotted. The hor- 
izontal axis is momentum transfer (Q) and the vertical 
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NiO, T=293 K, MARI, E =750 meV 
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FIG. 2: (Color online) Intensity map of expanded energy and 
momentum space. The solid white line would correspond to 
hydrogen recoil. 



axis is energy transfer. We have subtracted the high-Q 
background from the raw data. It is believed that the 
high Q data mostly consist of contributions due to mul- 
tiple phonons, and the high-Q data were fitted to a form 
of A + BQ'^ and subtracted from the low-Q data as de- 
scribed in Ref. [l9|. One can clearly identify a strong 
spectral feature around the excitation energy of 1 eV, 
which is more or less dispersionless. The spectral weight 
on the other hand shows quite strong Q dependence and 
the intensity falls to the background level around Q=15 
A~^, which is clearly seen in the integrated intensity plot 
shown in Fig. 1(c). The intensity at higher momentum 
is due to the background fluctuation. In addition to the 
1 eV feature, one can also notice a faint excitation fea- 
ture around 1.6 eV. This can be more clearly seen in Fig. 
1(d), in which we plot the intensity integrated up to Q= 
20 A-i as a function of energy. The well-defined peak 
centered around 1 eV and additional intensity around 1.6 
eV are identified as d-d excitations. 

To check if there are hydrogens in the sample, we mea- 
sured potential hydrogen recoil line [9]. In Fig. 2, we 
plot the energy and momentum region where hydrogen 
recoil line is expected (white solid line). As shown in this 
figure, there is no signature of any hydrogen scattering, 
eliminating the possibility that the observed excitations 
come from hydrogen. 

In order to understand the scattering mechanism and 
the observed momentum dependence, we have carried out 
a cluster calculation. We consider the electronic system 
to be described by a local model ^] in which the rele- 
vant degrees of freedom are the d-electrons of the Ni^+ 
ion in an octahedral ligand field. Our choice of model 
parameters Hi [HI are consistent with previous theory 
and experiments on d-d excitations of NiO, and were not 
adjusted to fit our data [32|. The model includes ex- 
plicit electron-electron interactions and is solved using 
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the exact-diagonalization method to obtain ah the eigen- 
states and eigenenergies. 

The eigenstates and energies are used via Fermi's 
golden rule to obtain the double differential scattering 
cross-section for unpolarized neutrons as 

f 

where [Hj 

V = X (^^1 J2 ^'"^'^^ (Pi + X Q) l^o) , (2) 

Q = — kj is the momentum transferred to the elec- 
trons, uj is the energy lost by the neutron, and and 
En are the states and energies of the electronic system. 

The matrix elements in Eq. (|2]) have a significant de- 
pendence on the CF and SO coupling, which can be illus- 
trated by considering the simple limit where Q ^ (i.e., 
the magnetic-dipole approximation). In this approxima- 
tion the allowed excitations are connected to the ground 
state by the operator L + 2S, where L is the total or- 
bital angular momentum and S is the total spin. For an 
isolated non-relativistic atom neither L nor S can cause 
inelastic transitions. However, in the presence of a CF 
it is well known that the diagonal matrix elements of L 
are "quenched" . The weight that is removed from these 
diagonal matrix elements effectively appears in the off- 
diagonal matrix elements; the presence of a CF allows L 
to contribute to the inelastic scattering signal. Similarly, 
the presence of SO coupling allows S to contribute to 
the inelastic signal. The magnetic dipole approximation 
illustrates the importance of CF and SO. 

However, due to neutron kinematics, the momentum 
transfers at energy-loss appropriate for d-d excitation (on 
the order of eV) can be large compared to the inverse 
radial size of the electronic orbit als. Thus the magnetic- 
dipole approximation is not always valid and the full ex- 
pression, Eq. (j2j), must be retained. We will see that the 
dipolar terms give rise to a non-dispersive intra-multiplet 
peak at ^ 10 Dq which dominates the spectrum at low 
momentum transfer. For higher momentum transfers, 
non-dipole excitations appear at higher energies (e.g., 
1.6 eV) and eventually dominate the relative scattering 
intensity. The origin of these peaks can be understood 
roughly from the Tanabe-Sugano diagram [22] for a 
system, which describes the change in multiplet ener- 
gies as a function of CF strength; transitions from the 
Ai octahedral symmetry ground state to the higher en- 
ergy T2 and Ti character states give rise to the spectral 
peaks at ~ 1.0 eV, and ~ 1.6 eV, respectively. The lo- 
cations of these peaks are slightly shifted and split by 
SO coupling. The relative intensity of each peak changes 
with the magnitude of Q as determined by the matrix- 
elements of Eq. ([2]). 
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FIG. 3: (Color online) Panel a) illustrates a false contour 
plot taken on MAPS with the Ei=500 meV. b) the intensity 
of the peak at 96 meV as a function of momentum transfer. 
The solid and open points are from MARI (Ei=750 meV) 
and MAPS (Ei=500 meV) respectively. The solid red line is 
the free ion Ni^^ form factor, c) The momentum integrated 
scattering and solid line is the calculated density of states 
from Ref. 



In Fig. 1(b) we show the theoretical INS intensity map. 
The energy integral from 0.9 to 1.25 eV over this map, 
and the momentum integral for Q < 20 A~^, are also 
shown in Fig. 1(c)- (d) as solid lines. The agreement be- 
tween theory and experiment is quite good, given that the 
only adjustable parameter is the overall amplitude. In 
particular, the momentum dependence shown in Fig. 1(c) 
suggests that the 1 eV transition occurs via magnetic 
dipole operator, L + 2S. The 1.6 eV transition spectral 
weight is only non-zero in the intermediate Q range, and 
vanishes both at small and large Q. Since SO coupling 
in this compound is small, the multiplet splitting due to 
the SO coupling only shows up as broadening of the 1.6 
eV peak. 

In addition to the high energy spectrum, we also 
measured the low energy spin excitation spectrum with 
Ei=500 meV, as shown in Fig. 3. A clear magnon mode 
is observed around 100 meV, which is the zone bound- 
ary spin wave energy. The assignment of this feature is 
unambiguous in INS. Earlier neutron scattering studies 
of spin wave dispersion reported that the zone boundary 
magnon energy is 117 meV [23|, [l^. In Fig. 3(b), we 
plot the intensity as a function of Q. The solid line is the 
magnetic form factor from the spin angular momentum 
contribution as tabulated in Ref. In Fig. 3(c), the 
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low-Q intensity (integrated up to 5 A~^.) is plotted as a 
function of energy transfer. Since we are using a powder 
sample, this integration is a quick method to obtain the 
magnon density of states. We compare this result with 
the calculated rnagnon density of states from Hutchings 
and Samuelsen [23|, which agrees very well. 

These lower-energy results shed light on understanding 
excitations observed with different spectroscopic tech- 
niques in the energy range of 100-200 meV. As shown 
in Fig. 3, no significant scattering intensity is observed 
above the one magnon band up to 250 meV. This should 
be compared with earlier photon spectroscopy investiga- 
tion. In optical spectroscopy studies [26|, a two- magnon 
mode at 193 meV was observed with Raman scattering, 
while an infrared absorption mode due to two-magnon 
plus phonon was observed at 250 meV. In a recent Ni 
Ls-edge RIXS study the excitations at 95 meV and 
190 meV were associated with transitions between atomic 
levels split by an exchange field . We note that "two- 
magnon" excitations in neutron scattering usually refers 
to longitudinal spin fiuct nations, which usually forms a 
broad continuum and has very small intensity [28]. The 
apparent lack of excitations around 190 meV in our INS 
data clearly illustrates the difference between the "two- 
magnon" excitations detected with photon spectroscopy 
and neutron spectroscopy. Further theoretical calcula- 
tion is necessary to address this issue quantitatively [29| . 

In summary, we observed d-d excitations of 1 eV and 
1.6 eV in NiO using inelastic neutron scattering. Such 
local multiplet transitions are allowed in the magnetic 
dipole and higher order channel due to the presence of 
non-zero crystal field and spin-orbit coupling. The ex- 
citation energy matches well with previous experiments, 
and the momentum dependence of the spectral weight 
could be well described by the cluster model calculation. 
Our observation illustrates that inelastic neutron scat- 
tering at high energies is a reliable and useful technique 
for probing electronic excitations in materials with sig- 
nificant crystal field and spin-orbit coupling. 
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